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The anisotropic superconducting properties of a Rb ;E Fe2-!/Se2 single crystal with T c ~ 32 K 
were investigated by means of SQUID and torque magnetometry, probing the effective magnetic 
penetration depth A e ff and the penetration depth anisotropy 7a. Interestingly, 7a is found to be 
temperature independent in the superconducting state, but strongly field dependent: 7a (0.2 T) < 4 
and 7a (1.4 T) > 8. This unusual anisotropic behavior, together with a large zero-temperature 
A c ff (0) ~ 1.8 /im, is possibly related to a superconducting state heavily biased by the coexisting 
antiferromagnetic phase. 

PACS numbers: 74.70.Xa, 74.25. Ha, 74.25.Bt, 74.25. Op 



I. INTRODUCTION 



With the discovery of superconductivity in 
LaFeAsOi-^Fj,, 1 a new family of iron-based high- 
temperature superconductors was found. Its simplest 
member is FeSei-^, which consists of a stack of FeSe 
layersP Its superconducting transition temperature 
T c ~ 8 K increases drastically with external pressure, 
reaching T c (8 GPa) ~ 36 Kj=™ Interestingly, a similar 
high T c ~ 30 K is attained in the iron-selenide family 
A x Fe2- y Se2 by intercalating alkali atoms (A = K, 
Rb, Cs) between the FeSe layers However, T c is 
found to decrease with pressure and is fully suppressed 
at 9 GPa for K x Fe 2 - y Se 2 (Ref. and at 8 GPa 
for Cs x Fe2_j / Se2P The critical temperature is almost 
insensitive to pressure below 1 GPa, suggesting that 
T c is almost independent of small variations of the 
lattice constants. This provides an opportunity to 
study the temperature dependence of physical quan- 
tities without being affected by changes of the lattice 
constants due to thermal expansion. Early /iSR experi- 
ments on Cs a; Fe2_ySe2 indicate that superconductivity 
and antiferromagnetism coexist on microscopic length 
scales.® The Neel temperature T N ss 500 K (Ref. HUJ) of 
Ko.8Fe2-j,Se2 is substantially higher than T c ~ 30 K. 
Several experiments point towards nanoscale phase 
separation between superconducting, vacancy-disordered 
domains and vacancy ordered antiferromagnetic (AFM) 
domainsP^Sl In contrast to the slightly hole doped 
FeSei^a,, the intercalation of alkali ions A into the 
FeSe structure introduces a large amount of electrons 
into the system.^ Other experiments^ suggested that 
this highly electron doped system contains no hole-like 
sheets at the Fermi surface, and thus electron scattering 
between hole and electron-like bands is impossible. 1 ^ 
Moreover, Fe vacancies in the crystalline structure 
order at ~ 600 KP^ In K ;E Fe2-, / Se2, the presence of 
vacancies appears detrimental to superconductivity.^ 



This intriguing microscopic coexistence of vacancy 
ordering, antiferromagnetism, and superconductivity in 
Aa;Fe2-j,Se2 points to an unconventional thermodynamic 
behavior of the superconducting state. 

Recently, the lower critical field H c \ was investigated 
in tetragonal K :E Fe2-j,Se2 for magnetic fields H along 
different crystallographic directions, i .e. pa rallel to the 
afr-plane and parallel to the c-axisf^EHl revealing a 
surprisingly low and isotropic (j,qH c i ~ 0.3 mT. Invoking 
the phenomenological relation between H c \ and the 
effective magnetic penetration depth A e g- for an isotropic 
superconductor fi a H cl = $ Q (lnK + 0.5)/(47rAg ff )p3 
a value of A e ff(0) ~ 1.6 — 1.8 ^m is obtained (as- 
suming an approximate Ginzburg-Landau parameter 
k ~ 100 - 200). This low field estimate of A cff de- 
viates remarkably from the small in-plane magnetic 
penetration depth A a &(0) ~ 0.29 [im derived from NMR 
experiments at 8.3 T on a similar sample.^ Moreover, 
the numerous observations of an anisotropic vortex 
lattice in high magnetic field d 19 * 2 2 ^^ contrast with the 
isotropic behavior of H c \ in low magnetic fields. In 
order to illuminate this intriguing field dependence of 
the superconducting properties, we performed a detailed 
magnetic study of Rb a: Fe2-i / Se2. 

II. CRYSTAL GROWTH 

A Rb a; Fe2- !/ Se2 single crystal with composition 
Rbo.77(2)Fei.6i(3)Se2 as refined by X-ray fluorescence was 
grown from a melt by the Bridgman method, using a 
presynthesised ceramic precursor of FeSeo.gs and metal- 
lic rubidium. For the precursor synthesis high purity (at 
least 99.99% Alfa) powders of iron and selenium were 
mixed in the molar proportion 1 Fe : 0.98 Se and pressed 
into a rod. This nominal stoichiometry, chosen on the ba- 
sis of our previous studies^ of the Fe-Se chemical phase 
diagram, provides an iron selenide of pure tetragonal 
phase. The mixture was prereacted in a sealed silica am- 
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FIG. 1. (Color online) Zero-field cooled measurements of 
m/H of single crystal Rb a; Fe2-j ; Se2 for (a) _ff||c-axis and (b) 
H\ |afe-plane. The change of vortex penetration with temper- 
ature is very similar for both directions, consistent with an 
isotropic (J,oH c i < 0.3 mT. The sharp transition at T c ~ 32 K 
demonstrates the high quality of the crystal. Due to demag- 
netization the magnitude of m/H varies by a factor of ~ 3 for 
both orientations. A rough estimation of the demagnetization 
factors from the sample dimensions yields 0.05 < TV 11 " 6 < 0.3 
and iV l|c ~ 0.65. 



poule at 700 °C for 15 hours and then grounded in an 
inert atmosphere, pressed again into a rod, sealed in an 
evacuated double wall quartz ampoule and resintered at 
700 °C. After 48 hours the furnace was cooled down to 
400 °C and kept at this temperature for 36 hours more. 
For the single crystal growth a piece of the Fe-Se rod was 
sealed in an evacuated silica Bridgman ampoule with an 
appropriate amount of pure alkali metal placed in an ad- 
ditional thin silica tube; 5% excess of Rb was added to 
compensate its loss during synthesis. The Bridgman am- 
poule was sealed in another protecting evacuated quartz 
tube. The ampoule was heated at 1030 °C for 2 hours 
for homogenization, followed by cooling down the melt 
to 750 °C at a rate of 6°C/h. Finally, the furnace was 
cooled down to room temperature at a rate of 200°C/h. 
After synthesis the ampoules were transferred to a He- 
glove box and opened there to protect the crystal from 
oxidation in the air. 



III. MAGNETIC MEASUREMENTS 

The superconducting properties of the plate-like crys- 
tal of dimensions ~5xlx 0.2 mm 3 (thickness of 0.2 mm 
along the c-axis) were characterized with a Quantum De- 
sign MPMS XL SQUID magnetometer. The temperature 
dependence of m/H, where m is the magnetic moment 
is shown in Fig. [I] for various magnetic fields H applied 
after zero-field cooling. The onset transition tempera- 
ture for this sample is estimated to be T c ~ 32 K. For 
both studied orientations (parallel to the c-axis shown in 



Fig.[lji, and parallel to the ab- plane shown in Fig.[l]3), the 
magnetic properties are very similar. For fioH = 0.1 mT 
a sharp diamagnetic transition is observed for both field 
orientations. In higher fields, the diamagnetism is rapidly 
suppressed, indicating that /io-ff c i for Rba;Fe2-ySe2 is 
very low (/io-Hci ^5 0.3 mT, including demagnetizing field 
correctionPSJ) 

A small piece of of rectangular shape and approx- 
imate dimensions 150 x 150 x 90 ^m 3 was cleaved off 
the above Rb;rFe2-ySe2 single crystal. Magnetic torque 
investigations were carried out using a home-made mag- 
netic torque sensorP^ For a uniaxial superconductor like 
tetragonal Rb K Fe2-j,Se2, the angular dependence of the 
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FIG. 2. (Color online) Angular dependent torque of single- 
crystal Rb :E Fe2-ySe2 measured at 15 K and 1.4 T. (a) Com- 
parison of unshaked and shaked magnetic torque. The un- 
shaked red data r cw (blue data r ccw ) are gathered by turning 
H clockwise (counter-clockwise) around the sample. The inset 
explains schematically the field configuration during the ex- 
periment. The shaked data are obtained by applying a trans- 
verse AC field, yielding reduced irreversibility and enhanced 
data quality, (b) Angular dependence of the shaked magnetic 
torque averaged for both directions. The green curve is the 
superconducting component of the signal, obtained after sub- 
tracting the background as explained in the text. A fit by 
Eq. |l} yields 7 A (15 K, 1.4 T) ~ 8.6 (black line). For compar- 
ison, the red dotted line is calculated with a fixed 7a = 3. 
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magnetic torque f — rh x in the London approxi- 
mation (H c i < H < H C 2) can be written according tcP^ 

(1) 

Here 9 is the angle between H and the crystallographic c- 
axis, e(0) = [cos 2 (6*) + 7^ 2 sin 2 (0)] 1 / 2 is the angular scal- 

I |c 

ing function, i?^ is the c-axis upper critical field, and r/ is 
a dimensionless parameter of the order of unity. Without 
loss of generality n is restricted to 1 within this workP^l 
The anisotropy parameter j\ = X c /X a b is the ratio of the 
out-of-plane and in-plane magnetic penetration depths, 
whereas A ff = (A^Ac) 1 / 3 denotes the effective magnetic 
penetration depth. It is possible that our sample presents 
phase separation between superconducting non-magnetic 
regions and antiferromagnetic, non-superconducting re- 
gions, with a domain size smaller than the penetration 
depth.EQ In that case, applying the Kogan model yields a 
A ff which may be renormalized to a larger value than the 
superconducting parameter (A^Ac) 1 / 3 . The field pene- 
trates more easily in this phase separated material, and 
the Kogan model yields an "averaged" effective bulk pen- 
etration depth. The good agreement between Eq. ([!]) and 
the data (Fig. [2]b) shows that the Kogan model is still 
useful, albeit with a broader interpretation of its param- 
eters. 

The magnetic torque experiments at various T and H 
were performed by turning H around the sample in a 
plane containing the c-axis (see inset of Fig. [2^,), and mea- 
suring the resulting torque. As an example, Fig.[2^L shows 
the angular dependence of the torque signal measured at 
15 K in 1.4 T. Note that the torque signal is affected 
by an angular irreversibility between the clockwise (r cw ) 
and counter-clockwise (r ccw ) branches. Such irreversible 
angular dependent torque signals are usuall y obse rved in 
hard superconductors due to vortex pinning!- 30 * 31 ! Due to 
the tetragonal structure of Rb 2: Fe2-ySe2, twinning in the 
crystal manifests itself in the results only through pinning 
on the twin boundaries, and all in-plane parameters are 
not separated along a- and b- axes in the analysis. In 
this work, the so-called vortex-shaking technique^ was 
successfully applied to reduce irreversibility, allowing a 
more reliable determination of the superconducting pa- 
rameters. This was done by applying a small AC field 
orthogonal to H in order to enhance the vortex relax- 
ation towards thermodynamic equilibrium. As seen in 
Fig. [2^, the vortex shaking clearly reduces the irreversible 
component, especially for H close to the a6-plane. In 
Fig. the average torque r = (r cw + r ccw )/2 is pre- 
sented. The torque signal consists of a superconducting 
component described by Eq. ([I]) and a magnetic back- 
ground component r BG = ~(Xab ~ Xc)V[i H 2 sin(20)/2. 
The variables \ c and Xab denote the magnetic suscepti- 
bilities along the crystallographic axes. Here Xab > Xc 
(as also mentioned in Ref.l22|) and tbg is large, consistent 
with a bulk antiferromagnetic phas e hav ing the magnetic 
moments aligned along the c-axis The subtraction 



of this antiferromagnetic background can be either per- 
formed by adding a sinusoidal component in the fitting 
routine of Eq. ([T]) , or by directly removing the symmetric 
sinusoidal component of the data as discussed in Ref. |3~T1 
All results presented in this work are independent of the 
background treatment. 

lie 

The parameters H]} 2 , j\, and A e ff can be extracted si- 
multaneously by analyzing the magnetic torque data with 
Eq. ([!]). However, in order to reduce the amount of free 

fit parameters, H}}£ was fixed according to a Werthamer- 
Helfand-Hohenberg (WHH) temperature dependence.^ 

The slope dH^/dT at T c was fixed to 1 T/K, in 

concordance with resistivity results of a similar sample.^ 
Small variations in dH^/dT do not affect the results 
of the analysis, since H C 2 contributes only logarithmi- 
cally in Eq. ([!]) and has no weight in the determination 
of 7a and A e g in low magnetic fields.^ Magnetic torque 
curves with the antiferromagnetic background subtracted 
are presented in Fig. [3] for various T (Fig. [3^) and H 
(Fig. [3)3). The insets show the normalized magnetic 
torque r norm = t(&)/ max[r(0)] close to the a&-plane. A 
change of the shape of T norm (6) qualitatively reflects a 
change of 7a . Note that the slope of r norm vs. changes 
strongly with H, but not with T. This demonstrates that 
7a is field dependent, but not temperature dependent. As 
the field direction is approaching the a6-plane, the screen- 
ing currents start to flow not only in plane, but also out 
of plane, which makes the torque depend strongly on X a b 
and A c in this angular region. The temperature and field 
effect on the anisotropy 7a are consequently already visi- 
ble on the data taken around the afr-plane, independently 
of the strict validity of Eq. in a phase separated ma- 
terial. 

The temperature dependence of A e g presented in 
Fig. [4^, can be fitted with the empirical power law 
A~ ff 2 (T) = A e - ff 2 (0)[l - (T/T c )"] with T c ~ 32 K, n ~ 5.2 
and A c ff(0) ~ 1.8 fim. This rather large value of the 
effective penetration depth may be due to phase separa- 
tion, as mentioned above. The value of n is quite large 
compared to what is expected from the two-fluid model 
the fitting expression comes from. This might be related 
to phase separation, as the links between superconduct- 
ing domains could depend on temperature, changing the 
temperature dependence of the penetration depth. Fig- 
ure [4Jd shows 7a (T) for various H. It is almost tem- 
perature independent between T c and T c /2 for all fields 
studied. The slight drop observed at higher tempera- 
tures in 1.4 T may be due to the proximity of the tran- 
sition, which decreases the superconducting signal. The 
increase of 7 at 1.4 T at low temperatures is most likely 
due to pinning effects.^ Most importantly, a remarkable 
monotonous field dependence of 7a is observed (Fig. |4jl) , 
with the strongest dependence in the lowest fields. In low 
fields, 7a(0.2 T) ~ 3.5. An extrapolation of the measured 
values of 7a towards zero field yields 7a (0 T) ~ 1 — 2. In 
Fig. [1J: the field dependence of A e ff(-ff) is shown for var- 
ious T. The effective penetration depth X c g(H) tends to 
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FIG. 3. (Color online) Angular dependence of the supercon- 
ducting component of the magnetic torque of Rb I Fe2- !/ Se2. 
(a) Magnetic torque at 17 K for various magnetic fields, (b) 
Magnetic torque at 1.4 T for various temperatures. The in- 
sets in both panels show the evolution of r norm close to the 
a&-plane. 



a constant value for /iqH > 1 T for all studied T. 



IV. DISCUSSION 

The extrapolated A c fr(0) ~ 1.8 /xm is surprisingly large 
compared to that of other iron-based superconductors. 
In the related iron selenide FeSeo.sTeo.s a much lower 
value of A e ff(0) ~ 0.7 /xm was reported.^ However, 
for Rb :E Fc2- a Se2 such a large A o g(0) is consistent 
with the very small noH c i J$ 0.3 mT observed in 
this work and in K :r Fe2_, / Se2 by others.^ The small 
7a (0 T) ~ 1 — 2 at very low fields is also in agreement 
with an isotropic H&. The increase of 7a = A c /A a b 
with H and the field independent A ff = (A^Ac) 1 / 3 

— 1/3 

imply that \ a b — A c fr7A decreases with increasing 
H. This is consistent with the high-field NMR result.^ 
The small high-field value of X a b reported in Ref. [22j 
in combination with the saturating A e g(J?) and the 
field dependence of j\(H) observed here, suggests that 
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FIG. 4. (Color online) Summary of all the results obtained 
by analyzing the experimental torque signal of Rb x Fe2-^Se2 
with Eq. M. (a) Temperature dependence of A^. The line 
is a power law fit to the data at 1.4 T. (b) Temperature de- 
pendence of 7a for various fields, showing that 7a is strongly 
increasing with H, but is almost independent of T. The dot- 
ted lines represent the average 7a for each field, (c) Field 
dependence of for various temperatures, (d) Field de- 
pendence of 7a for various T. The black line is a guide to the 
eye. 



7a continues to increase at higher fields. Note that 
the anisotropy of ~ 3 (Ref. |23|) observed in very high 
fields by upper critical field measurements cannot be 
directly compared with the magnetic penetration depth 
anisotropy 7a investigated here, because in general 

A field dependent 7a might be associated with a 
complex band structure, since in the case of multiple 
superconducting gaps originating from different bands 
the superconducting screening currents, related to A a ;, 
and A c , may give rise to an unusual behavior of 7a. 
A similar behavior was observed in MgB-2, where the 
two-gap excitation spect rum yields a strongly field de- 
pendent anisotropy ! 36 * 37 ^ In MgB 2 , as the field increases, 
the gap from the 3D cr-band is closed, so the main part 
of the superconducting fluid density comes from the 2D 
7r-band. The 2D character of the remaining band implies 
a larger anisotropy. 

Band structure calculations^®' in A x Fe2- y Se2 yield 
multiple bands, with four 2D sheets on the sides of the 
Brillouin zone and one cylindrical sheet at the center. 
Depending on doping, this cylinder can be split into 
two 3D cones in Ko.8Fe2Se2j^ which are completely 
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detached (3D) in the calculation of Ref. although 
when the authors use experimental lattice parameters 
these cones are replaced by a cylinder (2D). According to 
Ref. 201 the center band is a cylinder. In CsFe2Se2 and 
Cso.8Fe2Se2, this band consists in two almost detached 
conesP^In Tl .58Rbo.42Fei.72Se2p-^ two gaps of different 
amplitudes were observed on the side and center bands, 
but the candidate for a 3D inner band was to small for a 
gap to be observed. It is therefore not yet clear whether 
for Rb^Fe2-j,Se2 with the doping used in this work 
a 3D band is present. Clearly, more material specific 
work is needed in order to clarify the interplay between 
multiband superconductivity and the anisotropy of 
Rb a; Fe2- ! ySe2. However, if a 2D-3D band scenario is 
the origin of the field dependence of the anisotropy of 
Rba;Fe2-ySe2, it must be temperature dependent as 
wellp^ which is not the case according to Fig. [iji. 

The temperature independent j\ suggests that the 
origin of its field dependence is not related to the 
superconducting gap energy, which is strongly temper- 
ature dependent in the examined temperature range. 
However, there is one energy scale which is almost 
temperature independent in the superconducting state: 
the Neel temperature Tn ~ 500 KpSI A superconductor 
coexisting with an antiferromagnetic phase is expected 
to behave in a peculiar way, although the temperature 
range studied here is far too low to excite any change 
in this strongly coupled antifcrromagnet. However, 
magnetic field modifications can lead to changes in the 
domain structure, and therefore changes of coupling 
between superconducting areas. This could result in 
variations of the "averaged" effective anisotropy. As 
was shown for various iron-based superconductors, 
the lattice parameters, in particular the pnictogen 
height in the iron-pnictides, are directly related to 
superconductivity.^ Importantl y, s uch scaling also 
works for the iron-selenide layerPSl It is also possible 



that magnetostrictive effects, which are expected to 
increase with magnetic field, may influence the lattice 
parameters and by that the anisotropy of the system. 
In such a case the strong curvature of H C 2(T) in the 
vicinity of T c often observed in pnictides/chalcogenides 
may be explained by the change of T c with H caused by 
the change of pnictogen height. 



V. SUMMARY 

In summary, we present an investigation of the 
magnetic properties of Rb 2; Fe2_j ; Se2 revealing a strong 
field dependence of 7 a ranging from 7a (0.2 T) < 4 to 
7a (1.4 T) > 8. This behavior stands out among other 
iron-based superconductors, consistent with the singular 
coexistence of magnetic ordering and superconductivity. 
In accordance with lower critical field measurements, 
our data suggest that in very low fields 7a (0 T) ~ 1 — 2. 
At 1.4 T the effective magnetic penetration depth is 
A c ff(0) ~ 1.8 /im. The vortex phase in Rb 2 ,Fe2-2 / Se2 
is best described by an isotropic three dimensional 
state in low fields which becomes strongly anisotropic 
with increasing field. In this respect the novel iron- 
selenide Rb x Fe2-ySe2 could be a potential candidate for 
magnetic-field tuned applications of superconductivity. 
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